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Abstract. We study different one-parametric models of type la supernova magnitude evolution on cosmic time scales. 
Constraints on cosmological and supernova evolution parameters are obtained by combined fits on the actual data coming 
from supernovae, the cosmic microwave background, and baryonic acoustic oscillations. We find that the best-fit values imply 
supernova magnitude evolution such that high-redshift supernovae appear some percent brighter than would be expected in a 
standard cosmos with a dark energy component. However, the errors on the evolution parameters are of the same order, and 
data are consistent with nonevolving magnitudes at the lcr level, except for special cases. 

We simulate a future data scenario where SN magnitude evolution is allowed for, and neglect the possibility of such an evolution 
in the fit. We find the fiducial models for which the wrong model assumption of nonevolving SN magnitude is not detectable, 
and for which biases on the fitted cosmological parameters are introduced at the same time. Of the cosmological parameters, 
the overall mass density £1m has the strongest chances to be biased due to the wrong model assumption. Whereas early-epoch 
models with a magnitude offset Am ~ z 2 show up to be not too dangerous when neglected in the fitting procedure, late epoch 
models with Am ~ s/z have high chances of undetectably biasing the fit results. 

Key words. Cosmology: Cosmological Parameters, Cosmology: Observations, Stars: Supernovae: General, Surveys 



1. Introduction 

An important cosmological probe and the historically first 
indicator of a presently accelerating exp ansion of the uni - 
verse are type la super novae (SNe la, see Rie ss et al.l 4X998); 
IPerlmutter et alj dl999b ). SNe la are rare, massively luminous 
cosmic events that are believed to be violent thermonuclear ex- 
plosions of degenerate white dwarfs, resulting from the ignition 
of de generate nuclear fuel in stellar material dHovle & Fowled 
H%oh . Observations of nearby extragalactic SNe la indicate 
that they could be standardized, and could therefore be used 
as a distance measure. In a universe described by a Robertson- 
Walker metric, 

ds 2 = -c 2 dr 2 + cl 2 (t) [d X 2 + f(x)(de 2 + (sin 6) 2 d(f> 2 )] , (1) 

cosmic events of intrinsic luminosity _£ would yield a measured 
flux, 



on Earth, where ciq = a (to) is the cosmological scale fac- 
tor at the moment of reception (0), x - Xe ~ Xo the comov- 
ing distance between the emittor (e) and the receptor, and 
f(x) = {shut, x, sinhjt) according to a {negative, zero, pos- 
itive) curvature index k of the 3-space, respectively. Further, 
Z = (/le - ^o) Me is the redshift of the emitted light as observed 
on Earth. Equation (O accounts for the diminishing luminos- 
ity due to universal expansion of the emitted light-sphere, and 
to doppler-shift, where each contributes a factor (1 + z), and 
(1 + z) aof (x) is the luminosity distance that can be expressed 
as a function of redshift z: 



, c (1 + z) 

adz) = — r— f 




(3) 



r = 



(2) 



4tt(1 + z) 2 a Q z f(x) 

* Centre de Physique Theorique is UMR 6207 - 'Unite Mixte de 
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FRUMAM (FR2291). 



where Qk = —ktr/ (ciqHq) . Here, Hq is the Hubble constant. 
Via the Friedmann equations of General Relativity the Hubble 
function H (z) = HqE (z) can be expressed in terms of the dy- 
namical variables of the universe's matter and energy content: 



£ 2 (z;Q,-,w x ) 



Q k (l + zf + n M (l + z) 3 + Q R (1 + z) 4 + 



+ fi x Exp 



■I 



ln(l+ z ) 



(l+vv x (z))dln(l+z) 



(4) 
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Therein Q, depicts the matter densities of baryonic and dark 
matter (Qm = Qb + £2 C dm)> radiation (Qr), and the dark energy 
component (£2 X ), and w x = p x /(p x c 2 ) is the possibly time vary- 
ing equation of state of the dark energy component. The cur- 
vature parameter is determined by = 1-2 £2, . It is 
common to introduce the 'Hubble free' luminosity distance 
Dl = Hqc~ 1 cIl and to give (c/Hq) in Mpc. Adopting the defi- 
nition of apparent magnitude, m = -2.5 \ogT, we obtain from 
Eq.© the magnitude-redshift relation: 



m = -2.5 log T = M s + 5 log D L (z; Q,, w x ). 



(5) 



The normalization parameter M s = 2.5 log ^4n (c/Hq) 2 X" 1 ] is 
a constant number for all sets of luminous events with the same 
intrinsic magnitude M ma = -2.5 log X., and Eq.© becomes a 
fundamental and powerful tool for extracting the cosmological 
parameters {p;} = {£!;, vv x ) via the observables flux andredshift. 

The use of the simple relation Eq.© is complicated by sev- 
eral effects. The lightcurve data are not available over the whole 
frequency range, but only in a small bandpass that has to be re- 
lated to the whole bolometric luminosity. The adjustment to 
correct for this 'bandpass mapping' is done by adding a correc- 
tion term, whose value depe nds not only on t he source's red- 
shift but also on its spectrum ( Drell et al. 2000). The main com- 
plication, however, seems to be the deduction of the object's 
intrinsic brightness, which cannot be measured directly, from 
its spectroscopic and/or photometric properties. While SNe la 
clearly do not always have the same brightness, but indeed have 
a significant dispersion of 0.5mag in peak magnitude in most 
pessimistic estimations, one can establish a relation, however, 
between the decline rate Am^ (the total decline in brightness 
from peak to 15 days afterward) in the SN lightcurve shape and 
its absolute peak-brightness M pea k: 



Amis oc M, 



peak 



(6) 



0.15 dKowalski et al 



withi n an error margin of around cr M . 
2008b . which means distance estimates precise to w 7%. 
Several techniques of standardization have been developed to 
implement this relation int o SN data analysis, in particular 
the strech-factor method (Perlmutte 



template-fitting method (iHamuy et al.lll996allbF 7 The possible 



etalJ 1 19991) and the 



introduction of systematic errors by the use of diff erent fitting 



techni ques have been ex haustively studied by iDrell et al 
(2000) and lGuv et al] d2007l) . 



The use of Eq.© would furthermore be complicated if the ap- 
parent luminosity m was systematically varying with redshift 
due to various astrophysical effects. 

Intrinsic effects: A first class of possible effects would affect 
the intrinsic properties of the SN, like the luminosity M- mtx 
and the light-curve shape. First indications of this kind of 
evolution of SN properties came from the measurement of 
shorter rise-times preceding maximum lum inosity in SN light 
curves at higher redshift dRiess et alJI 1999D . Various possible 



reasons for a redshift evolution of the intrinsic SN properties 
have been proposed. Because the younger universe is charac- 
terized by lower overall metallicity, different metallicity of SN 



progenitors at different cosmic epochs has been considered 
as possibly causing an overall evolutionary effect. Also the 
white dwarf's initial mass and different mass fractions of the 
composing elements at the moment of explosion, varying 
with cosmic time, are suspect of affecting the intrinsic SN 
properties. The effect of varying composition and metallicity 
of the SN progenitors h a s been studied num e rically . The 
results of iHoeflich et all dl998l) . iLentz et alj d200ol) , and 



Ropke & Hillebrandt (2004) however disagree, in their predic- 



tions of Mj ntr and light-curve shape dependence on progenitor 
composition, but also use different types of explosion mech- 
anisms (detonation/deflagration models). Other effects might 
come from direct environmental phenomena like circumstellar 
gas clouds or magnetic fields. The progenitor's rotation is 
taken into account by three-dimens ional numeri cal models, but 
the results rest unconclusive so far (Ropke 2005). Furthermore, 
rather poorly understood dependencie s of the overall energy 
output on the progenitor's age remain. Gallagher et al. (2005) 
point out that progenitor age might more likely be the source 
of variability in SN peak luminosities than is metallici ty. A 



study of SN samples by the SN high-z team dStrolger et al 



2004) shows that SN models that require a large number 
fraction of prompt explosions (progenitor age < 2 billion 
years) poorly reproduce the observed redshift distribution 
and are rejected at > 95% confidence. When looking at the 
young universe, at redshifts z > 1.5, one clearly expects a 
younger mean progenitor sample, and a study of variation 
in SN magnitude with r espect to the mean p rogenitor age in 
different host galaxies bv lHowell et al.l ( 12007b indeed suggests 
more luminous SNe at higher redshift. 

Another problem with the progenitor age arises: the two 
principal models of supernova la progenitor systems, i.e. 
the single-degenerate and the double-degenerate progenitor 
system, differ in explosion delay time (and environmental 
conditions). The uncertainty in identifying t he progenitor 
system s may mimick an evolutionary effect ( Riess & Livid 
2006). 



Subpopulations JHamuy et alJ dl996allbh suggest that SN 
magnitude depends on the host galaxy (brighter SNe in early- 
type elliptic galaxies than in late-type spiral galaxies), and 
proposed a classification of SNe la into two subclasses. The 
study of lHatano et a l. (2000) also indicates two subpopulations 
that may be identified as resulting from different explosion 
mechanisms, i.e. plain detonation, deflagration, or delayed 
detonation models. Progenitor population drift leads to differ- 
ent SN subpopulations dominating at different cosmological 
epochs and may tr anslate into an overall evolutionary effect 
( Bran ch et al ]l200ll). The subpopulati on bias has rec ently been 
studied by Sarkar et al. ( 2008a ) and Lindei ( 20091) . showing 
that subpopulations with different mean peak magnitudes may 
bias the parameter of dark energy equation of state. 



Extrinsic effects: Another class of possible effects leading to 
an overall redshift evolution of apparent SN magnitude m 
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and light-curve shape consists of extrinsic effects, which do 
not come from the intrinsic SN properties being altered, but 
from astrophysical effects between the SN and the observer. 
Gravitational lensing effects have been considered to possibly 
mimick an SN magnitude evolution on cosmic time scales , 



but have been shown to be negligible dSarkar et alj [2008b) 



Also gray intergalactic dust has bee n pro posed to mimick 
an evolutionary effect since Aguirre (1999), and has been 
shown to possibly bias the extraction of the cosmological 
parameters dCorasanitill2006t Menard et al. 20091) . Also effects 
of the interstellar medium in the host galaxy may translate into 
evolutionary effects, the former systematically varying from 
early-epoch elliptic galaxies to late-epoch spiral galaxies. 



Given this rather rudimentary understanding of SN magnitude 
evolution in redshift range say < z < 1.7, we find it reason- 
able to test the optimistic assumption that the calibration rela- 
tions obtained for 'local', low -redshift SNe be without mod- 
ification valid for SNe at arbitrarily high redshift. We do not 
aim to extract evolutionary parameters or rule out models of 
SN magnitude evolution, but to study the systematics that are 
introduced when evolution is neglected in the analysis. 
To this end we developed detection criteria and tested their per- 
formance when confronted with the wrong model assumption 
on the SN magnitude. In section [2] we introduce some phe- 
nomenological models to describe systematical shifts in SN 
magnitudes, discuss general aspects of SN magnitude evolu- 
tion, and give limits on the model parameters obtained from 
current observational data in section[3] In the main part, section 
|U we study the implications, an unaccounted for evolutionary 
effect on SN magnitudes would have on data analysis. Therein 
we make use of a simulated SN sample consisting of ~ 2000 
SNe up to redshift z — 1.7. We present our approach to the de- 
tectability of such an effect, and discuss the possible introduc- 
tion of nondetectable biases on the cosmological parameters. 
From merged detectability and bias results, we introduce the 
concept of the danger of the phenomenological models. This 
paper ends with section [5] where the results are summed up 
and discussed. 

2. Parameterizations 

The possible evolution of SN magnitude with redshift is al- 
lowed for by a simple modification of the magnitude-redshift 
relation, i.e. in Eq.© we allow for a magnitude altering 

Am(z) evo , 



M s + 5 log D L (z; Q„ W x ) + Am(z) e 



(7) 



supposed to originate in an overall time evolution of the SN 
peak luminosity. In the following we simply type Am ev ° instead 
of Am (z) evo . Several phenomenological models for an SN mag- 
nitude evolution have been used in the literature. Because we 
know little about the physics of a magnitude variance with cos- 
mic time, these parameterizations somehow are a 'matter of 
taste' and are differently motivate d accor ding to the purpose 
of the respective analysis. iLinderl (120061) studies a Am evo oc 



zl (1 + z ) model for reason that i t is bounded at high redshift 
z — > oo. Ferramach o et al. I d2009l) use a Am™ oc (t (z) — to) in 



real data analysis, where t is the cosmic time: 



r(z; Q,, Wi) 



oo 



1 



(1 + z)E(z; O,-, wd 



dz. 



(8) 



In the ACDM case, Eq.([8j, however, is well approximated by 
a Am evo ~ -\fz model in the matter-dominated epoch z > z eq , 
where z eq is the redshift where Qm (z) - fix (z). In the dark- 
energy dominated regime z < z e q Eq.® is equivalent to 
Am™ ~ log ( 1 + z ) . Thi s logarithmic model was first consid- 



ered by lDrell et al.1 (I2000I) . who motivated it as a parameteriza- 



tion of the intrinsic luminosity X: 



(9) 



where fl\ is a real number. This assumption on a possible red- 
shift dependence of the intrinsic SN property trivially translates 
into the parameterization 



Am™ =2.501 log(l +z), 



(10) 



which will be referred to as Model 1 throughout this study. We 
in some respects prefer this consideration of a varying X in- 
stead of modifying m, for reason of the latter being an artificial 
construct. Since we would expect a cosmological drift of SN 
magnitudes at very high redshift to come from astrophysical 
effects more or less directly related to cosmological expansion, 
and assuming that X. itself is more or less a direct function of 
these variables, we would expect the effect on £ to be natu- 
rally related to the scale factor of the universe, a (t), and thus to 
(1 + z). Being aware of the parameterization dependence of all 
analysis results, we complement our study by considering the 
general models 



Am evo = p lZ a 



(11) 



to study the contribution of powers in z to our results, a and /?2 
again are real numbers. For a — > in Eq.dTTI) the evolution- 
ary effect becomes highly degenerate with M s . In particular, 
we use the three submodels a = 0.5, 1, and 2, which will be 
referred to as Model 2.1, Model 2.2, and Model 2.3, respec- 
tively. These three power contributions are chosen in order to 
give more weight either to low-z or to high-z SNe in the sam- 
ple. Model 2.1 would model a magnitude evolution dominant at 
lower redshift (z ^ 1), meaning a late-epoch model. Model 2.3 
would give a magnitude evolution dominating at high redshift 
(z ^ 1) and will thus be considered as a early-epoch model. We 
note that Model 2.2 had alr eady been th e point of a qualita- 
tive discussion in th e work of lRiessI (120001) . and was studied by 
Nordin et all (2008). We also note that a Model 2 with a * 0.6 
describes the functional form of Eq.® correctly in the whole 
redshift range < z < 1.7, and at the same time is a good 
approximation to the oc zl (1 + z) parameterization. 

3. Real data 

3.1. Constraints on evolution 

We obtained current constraints on the /3 parameters by 
a combined fit on the values of the cosmic microwave 
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background (CMB) shift parameter R = 1 71 ± 0.03 (see 
Efstathiou & Bondlll999tlWang & Mukherieell2007t) . the bary- 
onic acoustic oscillation (B AO) parameter A - 0.469 ± 0.017 



(see lEisen stein et al. 2005), and the 307 SN magnitudes as 
compiled by iKowalski etail d2008l) . Instead of a full CMB 
and BAO analysis, R and A are being used to save computing 
time, cf. section [4] We impose flatness 2,-Q,- = 1 on all fits. 
To take a possibly time-varying value of the equation of state 
of the dark energy component into account, we used the 
Chevallier-Po larski-Linder (CPL, Ichevallier & PolarskilEooT 



Linderl 120031) parameterization W x (z) — wo + w a z/(l +z), 



where w a = -2jj^| ^ depicts the overall time variation of 
w x . The fit parameters thus are {p,} = {M s , Qm, wq, w d ,ft}. We 
neglected the radiation component fi« in the fit. 

The fit was performed with the cosmological analysis tool 
gfitQ The best fit is found by a minimization of the 



X 



[m - m fit (z; Pi )f V- l [m- m fit (z; Pi )] , 



(12) 



where m - (mj,R,A) is the vector of the SN magnitudes 
mi = {mi, ...,m-}ffj] and R and A. The vector m nt (z;pi) is 
the corresponding vector of values for the fitted parameters, 
and the covariance matrix V is the inverse of the error matrix 
whose elements are made up of the errors <x, on SN magnitudes 
and on R and A: V,-, -1 = cr, 2 . The errors on parameters pk are 
obtained by minimization over all other respective parameters 
{Pi)i*k for p^s in an environment around the best-fit value p^ x . 
One sigma errors ±lo~ Pk then correspond to^. + 1. 



Table 1. Best-fit values for models 1 and 2. 



a) 


Q.M 


P 


w x = const. 


X 1 /" 


i) 

Model 1 
Model 2.1 
Model 2.2 
Model 2.3 


27 +U U2 

' -0.02 

24 +0 05 

"■^-0.04 

29 +0 05 
24 +0 04 

" ,z - -0.03 

26 +0 02 

"■^"-0.02 


-0.22!°;! 
+0.08+° 2 ° 
-0 15 +0 15 


-0 95+ UU8 

"■''-'-0.08 
1 1 O+0.34 
_1 - 1S -0.34 

-0.86!» 24 
-1.16t 22 

-l.oot! 2 


0.72 
0.71 
0.71 
0.71 
0.71 


b) 


Q.M 


p 




X 2 /n 


ii) 

Model 1 
Model 2.1 
Model 2.2 
Model 2.3 


n 90+O.O2 
U - 2 ° -0.02 

24 +0 04 

u - ZH -0.04 

23 +0 05 
25 +0 03 

"•^ -0.03 

26 +002 

" ,z - -0.02 


-0 24 +0 - 23 

"'^-0.23 

-0 30 +031 

u '- 3u -0.30 

-o.i3t!i 

-0-08^ 


1 1 4+0.16 n SP+ u - 4y 

_i 40+0.38 1 on+0.40 

1,HO -0.31 * -Q.6Q 

- 1 65 +0 - 55 1 72 +0 ™ 

1,UJ -0.47 1 -'- i -1.03 

-136+ 024 1 19+051 

l,J -0.24 1,1 -0.91 
1 91 +0.18 r> Qi +0.48 
1 - ii -0.18 yJ - yl -1.09 


0.71 
0.71 
0.71 
0.71 
0.71 



The results for all four models are assembled in Table[T] where 
we show the results for the case w x = const., and the dy- 
namical case w x = wo + w^zj (1 + z), case a) and b), respec- 
tively. The table also shows the results when /3 = was im- 
posed on the fit, cases i) and ii), respectively. These agree with 



1 GFIT is a cosmological analysis tool developed by A.T. that al- 
lows studies of real or simulated data for various probes (SNe, CMB, 
BAO, weak lensing (WL), clusters, ...). Information on the code and 
its utilization can be obtained by contacting: tilquin@cppm.in2p3.fr. 



the results of IKowalski et ail d2008l) . iKomatsu et ail d2009l) . 
and Tilquin et al] ( in preparation!) . In the following we call the 
model with {fi M , w , w a } = {0.28, -1,0} the ACDM model. 
Concerning the x 1 per degrees of freedom n of the fit, all mod- 
els of Table [1] perform equally well, i.e. there is no significant 
tendency toward considerably higher or smaller x 1 f° r an y of 
the models. 

When we impose w x - const., we see a surprisingly high best- 
fit value of y8i = -0. 22+^32 f° r Model 1, which corresponds 
to a magnitude drift Am evo |i. 7 = -0.24+^ at redshift z = 1.7. 
The Om best-fit value is decreased by 0.03 compared to the 
standard j3 — fit, which gives Om = 0.27. All {/?,}, how- 
ever, are consistent with ACDM at the lcr level. Model 2.2' s 
evolutionary parameter has a lower negative best-fit value and 
is consistent with non-evolutionary SN magnitudes at the \<r- 
level, while w x in this model is closer to the ACDM value -1, 
and errors on the parameters f3 and w x are smaller. Model 2.3 
again shows smaller errors on the parameters /3 and w x , and the 
Q.M and w x parameters approaching the ACDM values 0.28 and 
-1, respectively. Its evolutionary parameter [}, however, shows 
up to be inconsistent with non-evolving SN magnitudes at lcr. 
Model 2.1 is a special case in this table, since it is the only one 
with a positive best-fit value for ft. But it is consistent with no 
evolution at all at lcr. We find that at a « 0.6 the sign of the 
best-fit /3 turns around, and that we fit ft s* for all a ^ 0.6. Due 
to correlation between ft and w x , the best fit value of w x also 
changes qualitatively with decreasing a, and we fit w x s* -1 
for all a § 0.6. 

When we allow the dark energy equation of state to vary 
with time, i.e. if we include w- d in the fit, the picture changes 
considerably. In Model 1 the wq best-fit value drops signifi- 
cantly below - 1 . Equally the w a value is high, and the fit be- 
comes inconsistent with a cosmological constant {wo,w a } = 
{-1,0} at the lcr-level. For Model 1 the w a value even is incon- 
sistent with its cosmological constant value at 2cr. For models 
2 we find the tendency that higher powers in redshift results 
in lower \ft\ best-fit values, and Qm best fit-value increasing to- 
wards Q.m = 0.28. wo equally moves towards its ACDM value 
-1, but remains inconsistent with it at the lcr level. The CPL 
parameter w- d is inconsistent with w a = also in models 2. 1 and 
2.2, but becomes consistent with w. d = in Model 2.3. We note 
for Model 2.1 that the introduction of the w d parameter pushed 
the a = 0.6 limit, where the best-fit value of ft became positive 
and wo > -1 in the w x = const, case, to lower redshift powers 
a ~ 0.1 due to the strong degeneracy between ft and w a . All 
models including Model 2.1 have ft < best-fit values. 
Figure Q] illustrates the best-fit values of Table Q] in a magni- 
tude redshift diagram; i.e. , we plot the magnitude drifts Am evo 
obtained from the best-fit ft values via Eqs. ( TTOb and ( fTTT i over 
redshift. The figure highlights that all Am evo > are excluded 
at lcr when w x + const, in the fit. The data are consistent 
with Am evo = in all the fits, with two exceptions. Model 2.3 
exludes the hypothesis of nonevolving magnitudes at lcr in the 
case where a constant dark energy equation of state w x = const, 
is assumed. And Model 1 does the same when w a is included 
in the fit. 
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a) 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
w x =const. 

















0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 1,6 



b) 




redshift z 



Fig. 1. Illustration of the Icr limits for models 1 and 2 ob- 
tained by an SNe+7?+A fit in a magnitude vs. redshift diagram, 
a) Dark energy equation of state w x = const, has been imposed 
in the fit. b) w x = wo + w a z/(l + z) is allowed to be dynamical. 



Table[2]shows all contours we obtain by a Model 1 fit, with the 
four different fit-parameter sets i), ii), a) and b), respectively, 
cf. Table Q] Contours in parameter spaces (pj,pk) are obtained 
by a ^-minimization over all other parameters {ptit^k, and 
{68%, 95%) confidence level contours correspond to contours 
X 2 min + {2.3,5.99}, respectively. The parameter pairs pj x in 
the headings of the columns depict the abscissae x ordinates 
of the contour plots. We see that introducing the evolutionary 
parameter/?! in the fit results in a decrease of the best fit Qm 
value, which as we saw above is valid for all models whenever 
y6 fit < 0. Concerning the present value of dark energy's equation 
of state wo, including the f3i parameter in the fit has the same 
effect as including w fl , i.e. it moves the wo best-fit value beyond 
the phantom barrier w x = - 1 . Whereas including either w a or 
/3\ keeps the fit result consistent with wq = — 1, including both 
of them at the same time makes w$ inconsistent with Wq = - 1 . 
This effect persists in all models, but is the strongest in the late- 
epoch model 2.1 and the logarithmic model 1. 



3.2. Hubble diagram 

The general tendency to prefer a negative (negative magni- 
tude drift Ara evo ) and thus brighter SNe at higher redshift can 
also be seen from the Hubble diagram of the SNe sample. We 
bin the sample in redshift bins of width 0.1. Figure|2]shows the 
relative magnitude deviation 



Amu 



= — V (m - mh 

JVbm ,'ebin 



(13) 



of the SNe sample from the best-fit magnitudes of case ii). 
Here, M,in depicts the number of SNe in the bin. Error bars 



1 



= — 7(- 

M. ; „ n-„ 



(14) 



come from th e erro rs cr m . on magnitudes as reported by 



Kowalski et al.1 (120081) . We included in gray the distribution of 
SNe in the plot, the right ordinate showing the number A^in of 
SNe in the bin. High-redshift supernovae, i.e. SNe at z ^ 1, 
seem to have a tendency to appear more luminous than the best 
fit. A linear least residual squares fit on this binned magnitude 
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Fig. 2. The relative magnitude deviation Amb; n , Eq.(fT3b, of the 
SN sample from the magnitudes we obtain from the case ii) 
best fit. 



deviation gives a negative slope of 

Am evo 

= -0.04 ± 0.06 



(15) 



magnitude per redshift, which corresponds to a magnitude 
drift Am evo |i. 7 = -0.07 ± 0.10 at redshift z = 1.7. Data thus 
seem to prefer more luminous SNe at higher redshift, but are 
consistent with nonevolving SN magnitudes at Icr. We do, 
however, see clearly from Fig. |2]that with the current data the 
statistics is poor at high redshift, with only 18 SNe at z > 1. 
More high-redshift SN magnitudes are indispensable to check 
for m agnitude evolution, as frequ ently pointed out since iRiess 
(|2000h . (e.g. IRiess & Livioll2006h . 
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From sects. 13. I land l3.2l we draw the conclusion that, while 
preferring negative magnitude drifts Am evo |i 7 < 0, data show 
no d efinite indication of SN magnitude evolution, in ag r eement 
with iKowalski et al.1 d2008l) . iFerramacho et all (120091), iLinder 
(booC and the res ults reported by iBronder et al.l (120071) and 



Sullivan etal. (2009). 



4. Simulations 

We now study the impact, an unaccounted for SN magnitude 
evolution would have on the values of the extracted cosmo- 
logical parameters in the framework of a combined analysis. 
Our database is a set of 2000 simulated SN magnitudes in red- 
shift range < z < 1.7, from which ~ 1000 at z > 1. The 
sample corresponds to what could be expec ted from a satel - 
lite mission with good control of systematics ( Kim et al.l20 04). 
We added 300 SNe in the nearby range z < 0.1 as expected 
by nearby SN la survey s like the Nearby Supernova Factory 
( Wood -Vasev et al . 2004). The sample was binned into redshift 
bins of width 0.1. The intrinsic error on each SN magnitude is 
assumed to be cr intl - = 0.15, and systematic errors cr sys = 0.02 
on magnitudes are included in the analysis. This systematic er- 
ror enters the covariance matrix by Vu = <x, 2 + (cr,-) 2 ys and 
Vij = rij(<Ti) S y S (a-j) sys , where r,j is the correlation coeffi- 



cient between bins, which however is set to ry = 0. In the 
simulation of the SN magnitudes we allowed for a redshift 
evolution according to Eqs. (TTOb and (fTTT >. Scan range in (3 is 
-0.5 < p < 0.5 for Model 1 and Model 2.1, -0.4 < f3 < 0.4 
for Model 2.2, -0.3 < (3 < 0.3 for Model 2.3, and the 
other respective fiducial cosmological parameters are fixed to 
{fi M , wo, wa} Fid = {0.30,-1,0)0 We checked the stability of 
our results to variations in these fiducial parameters and found 
no significant variance from our general conclusions. Fit pa- 
rameters are the standard parameters {p,} fit = {M s , Qm, wq, w a }. 
The magnitude evolution inherent in the simulated data was 
thus neglected in the fitting procedure; i.e. , we imposed [3 = 
on the fit. To economize computing time we did not perform a 
combination with full BAO and CMB data, but simulated the 
CMB shift parameter/? and the BAO reduced parameter A. We 
find this simplification justified, because our task is not to ex- 
tract real cosmological parameter values but to study the sys- 
tematics introduced by a possibly wrong assumption on one of 
the probes. We assumed an error of ±0.01 on R (which is the 
estimate for future PLANCK data) and an error of ±0.005 on 
A, to represent an experimental setup available in the near fu- 
ture. 

The numerical tool for our parameter study with simulated data 
is KOSMOSHOW0 which implements the ^-minimization 
procedure as described in section [3] The errors on the cosmo- 
logical parameters {/;>,■ } are estimated at the minimum by using 
the first-order error propagation technique: U = J.V.J T , where 
U is the error matrix on the cosmological parameters, and J the 



2 In cases where a confusion is possible, we added a superscript 
'Fid' to fiducial parameters and a superscript 'fit' to fitted parameters. 

3 KOSMOSHOW is a cosmological analy- 
sis tool developed by A.T., and is available on 
http://marwww.in2p3.fr/renoir/Kosmo-Pheno.php3 



Jacobian of the transformation. The fit of the simulated data 
gives the central values {/?,} and errors {cr p .} for the parameters, 
along with their correlations p a ;,, which in general are strong. 
We neglect the radiation component Qr and continue assuming 
spatial flatness in the following, 2, flj = 1. 

4.1. Illustration 

We observe that, even with this high statistics at redshifts z > 1 
and forecasted small errors on R and A, an unaccounted for 
evolution of SN magnitudes will possibly lead to misinter- 
preting the fit results. Suppose a fiducial cosmological model 
{Q M ,wo,w a ,y0i) Fid = {0.255,-1.20,0, -0.10}. This isacosmos 
with a dark energy component with a constant equation of state, 
w x = const. = wo, where SN magnitudes evolve with red- 
shift according to Eq.dTOl). This fiducial model is represented 
by the small circle in the (Qm, h>o) parameter space of Fig. [3] A 
fit on this cosmology's SN magnitudes combined with R and 
A gives a best fit {D. M , W , w a } fit = {0.274 ± 0.007,-1.12 ± 
0.09,0.23 ± 0.41} with a very low^ 2 /n = 0.2 when the pos- 
sibility of an SN magnitude evolution is not taken into ac- 
count in the fitting procedure. Here, n depicts the number of 
degrees of freedom of the fit. The best-fit value and the cor- 
responding 68% and 95% confidence level contours are plot- 
ted in Fig. [3] which highlights that the fit result would lead us 
to reject the true cosmology at more than the 95% confidence 
level. We point out that the best fit on this cosmology is consis- 
tent with the actual real data best fit on SNe+BAO+CMB+WL 



(Tilq uin et aLllin preparat ion), which pins down the parameter 
values to {Q M , w , w a } fu11 fit ~ {0.28,-1.1,0.4}. Even with a 
mid-term prospective SN dataset and combination with R and 
A with forecasted small errors, we would therefore risk misin- 
terpreting our best fit because of a wrong model assumption on 
one of the probes. We have been unable to find an illustration of 
a confusion between a fiducial ACDM model with an SN mag- 
nitude evolution, and the 'full fit' model given above. Indeed, 



we find that w a fit < w a Fld (= 0) whenever w. 



Fid 



< — 1, whatever 



sign and model of Am ev0 . 

To avoid this misinterpretation, one of course would like 
to have the means in hand to possibly detect a wrong model 
assumption on one of the probes. We introduce three differ- 
ent generic criteria that allow judging the performance of a fit. 
We then check, for which fiducial cosmologies involving an 
SN magnitude evolution, parameterized by /3, these criteria are 
strong enough to detect the wrong model assumption in the fit. 

4.2. Statistical goodness criterion 

As a first criterion for judging the goodness of the fit we use 
the ^- 2 -test. To evaluate the statistical go odness of the fit, w e 
choose the^ 2 ut of the so-called 'p-value' (lAmsler et al ] l2008l) : 



1 - 0.68 (0.95) : 



r°° x nl2 - l e- xl2 



(16) 



as statistically justified cut-off values, and consider all models 
with^ 2 > xl u t as ruled out at 68% (95%) confidence level, re- 
spectively. Here, n is the number of degrees of freedom of the 
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Fig. 3. Illustration of the possible risks introduced by the 
wrong assumption of nonevolving SN magnitudes. The cir- 
cle depicts the fiducial cosmology {Qm, wq, w fl ,y8i} Fld = 
{0.255,-1.20,0,-0.10}. The best fit with the standard set of 
cosmological parameters, thus neglecting SN magnitude evo- 
lution, is depicted by the full circle, see text. Black lines are the 
68% and 95% confidence-level contours. 



fit. The detectability limits on the fiducial parameter /? we ob- 
tain by virtue of this criterion are tabulated in Table [3] We fur- 
thermore test for the detectability of an evolutionary effect on 
the SN magnitudes by a consistency criterion which is genuine 
to SN data. 



4.3. Consistency Criterion 

In the literature it is a common approach to consider the terms 
Tow redshift' and 'high redshift' when discussing SN data, 
where the definition of the limiting redshift between the two 
ranges, however, is not uniform. Some authors refer to the scat- 
ter of the Am in the Hubble diagram Fig. [2] and divide the 
sample into high and low redshift subsamples with equal mean 
scatter, while others cut the sample by simply assuming that 
z > 0.5 or z > 1 be high redshift. We find it a good statistical 
means to cut an SN sample into high and low redshift sub- 
samples at a certain z cut such that the errors on the obtained fit 
parameters pk be equal for the two subsamples. That is to say, 
we impose cr J" = <x J 1 on the high (hr) and low (lr) redshift 
subsamples. We therewith ensure that the two subsamples have 
equal statistical weight with respect to the parameter pk, but 
will find different cut redshifts for different pt- For our simu- 
lated SN sample as outlined above, we find z cut = 1.20 (0.95) 
with respect to equal errors on pk = Dm (wo). These z cut are 
found by looking at the SNe only (and not at the combined set 
SNe+7?+A), because we want it to be a criterion intrinsic to 
SNe. 

Figure [4] illustrates the application of this criterion for a 
simulated Model 1 dataset with -0.5 < /? Fid < 0.5. All other 
fiducial parameters are fixed to {Q M , wo, w a } Fld = {0.30,-1,0} 
in the simulation, and the SN sample has been split in two at 
redshift z cut = 1 .20. The figure shows the wo obtained by sepa- 



rate fits on the SN lr -kR+A and the SN hr +R+A datasets, plotted 
over the varying /? Fld . We recall that the fitted parameters are 
{Pi} & = {Ms, Qm, wq, w a ) as throughout this section. 




x— z<1 .2 subsample SN ' 
o— z>1 .2 subsample SN hr 



tl T T T T T 



-0,5 -0,4 -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 0.4 0,5 

fiducial evolution parameter p 

Fig. 4. Illustration of the consistency criterion. We plot the 
w Q -values obtained for the low redshift z = 0...1.2 (black 
crosses) and the high redshift z = 1.2. ..1.7 (blue circles) sub- 
sample over the fiducial evolution parameter /? Fld of Model 1 . 
Error bars are the lcr errors cr„! r and cr„, . 



We find, as should be expected, consistent results around 



0, i.e. l. the wi" obtained by a fit on the high-redshift 



subsample lie well within the errors of the w () h obtained by a fit 
on the low-redshift subsample, and vice versa, and 2. both the 
w hl and the w h are in agreement with the fiducial w Fld = - 1 . 
But the wo clearly become inconsistent for larger |y6 Fld |. To give 
the exact inconsistency bounds we look at the consistency bias 
C, = \pf lr - p^\ and consider 



Ci 



< i 



(17) 



<o 2 ■ 



• (<) 2 



as the criterion of consistency between the low and high red- 
shift subsamples. Here, T is a factor corresponding to the 
68% (95 %) confidence level. Values for T are tabulated in the 
literature dAmsler et al . 2008) or can be obtained from integrat- 
ing 1 - 0.68 (0.95) = j T T f(t, n)dt, where / is the Student's 
function and n the number of degrees of freedom of the fit. 
For the specific example of Fig. [4] we find from Eq.dTTt: 
-0.17 (-0.40) < Pi < 0.13 (0.23) as the fiducial parame- 
ter range where an SNe+/?+A analysis gives consistent wo-fit- 
values on the high- and low-redshift SN subsample. All limits 
on all models obtained by this criterion are tabulated in Table[3] 
This table also compares the results for z cut = 1 .20 with those 
for z cut = 0.95. 

4.4. Physical criteria 

We impose 

(w + w a )<0 (18) 
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on the fit results. Equation (fl~8T > is implied by the condition that 
w x be smaller than zero at the time of decoupling, w x (zcmb) < 
0, to allow for a matt er-dominated early un iverse and thus 
structure formation (see Steinhard t et al 
We also impose 



Q M > 



(19) 



on the fit to allow for the presence of matter in the universe. 
It however turns out that in the context of a combined anal- 
ysis SNe +R +A neither Eq.dT8b nor Eq.(fl9ll contribute any 
constraint on the yS ; Fld in the range of study. We therefore do 
not discuss these criteria in the following, and omit them from 
Table El 

4.5. Results 
4.5.1. Detectability 

Putting together the criteria developed above, we obtain /3- 
ranges for all four fiducial evolutionary models where an SN 
magnitude evolution would pass the fitting procedure unde- 
tected by these statistical and physical detectability criteria. 
Table [3] shows the ensemble of the results. The consistency 
check is particularly powerful, and on the 68% confidence level 
yet more decisive than the^ 2 -test in all cases. At the 95% con- 
fidence level, the x 1 -check performs better for (3 < for all 
models. We note that the consistency check for Qm is most 
powerful when the sample is split in two with respect to cr„ , 
i.e. when z cut = 0.95. Inversely, when the cut is done with re- 
spect to <tq m , i.e. at z cut = 1.20, the most constraining test is 
the consistency check on wo- We underline the most constrain- 
ing detectablitity limits in Table [3] We can translate the under- 
lined 68% confidence level /Mimits from Table [3] into magni- 
tude shifts Am 6 ™!!.? at redshift z = 1.7 via Eq. ([TOj and ( fTTT i. 
and obtain the Table [4] 



Table 4. The Am e 



[1.7 



■intervall in which an SN magnitude 



evolution would not be detectable at the 68% confidence level. 



Model 



1 

2.1 
2.2 
2.3 



Non-Detectability Zone 



[-0.18,0.14] 
[-0.25,0.20 ] 
[-0.15,0.12] 
[-0.12,0.12] 



We highlight that, at the 68% confidence level, the most per- 
forming detection criterion is the consistency check. However, 
as noted, when going to higher confidence levels, one may find 
the^ 2 -check better suited to detecting the effect of SN magni- 
tude evolution. From these Am™ limits one concludes that the 
power models Eq.dTTI) are the higher constraint, the larger the 
exponent a is. The early-epoch model 2.3 is detectable when- 
ever the evolutionary magnitude shift |Am evo | 1.7I > 0.12 at red- 
shift 1.7. The late-epoch model 2.1 would pass undetected for 
all evolutionary effects |A;« evo |i 7 | < 0.20 even in an SN survey 



with redshifts up to 1.7 and combination with R and A. The 
detectability limits derived in this section also underline the 
tendency that a negative magnitude drift Am evo < is harder to 
detect than a positive one. We saw the preference for Am™ < 
from real data in section|3] cf. Figs. Q] and [2] 

4.5.2. Biases 

Whether this undetectability of the SN magnitude evolution is 
dangerous or not depends on the quality of reconstruction of the 
fit parameters {/?,}. For a fiducial cosmology involving, e.g. , a 
j3\ = -0.03 (corresponding to Am ev °i 7 = -0.03) evolutionary 
effect on SN magnitudes, the problem would not be detected. 
But we would not need to worry because all the fiducial pa- 
rameters {pi} Fld = {Qm, vvo, w a } Fld lie within the lcr error of 
the fitted value, and the fiducial cosmology would thus be 'cor- 
rectly' reconstructed at the lcr level despite the undetected SN 
magnitude evolution. We introduce the bias S,- on parameter pi 
as Si = [p, - pf I, and adopt the notion that the parameter pi 
is 'biased' whenever: 



* >1. 



(20) 



By applying this notion, in this section we restrict ourselves to 
the 68% confidence level. We call the zones of fiducial param- 
eter space where parameters p^ are reconstructed with biases 
the bias zones. Correspondingly, the zones of correct recon- 
struction without bias are referred to as validity zones. 
Whereas no bias is introduced in the fit in the preceding exam- 
ple, fiducial models exist where an evolutionary effect Am em 
would pass our criteria undetected, and where, at the same time, 
one or several of the {/?,} are biased. Our simulations, e.g. , 
show that in a universe with an unaccounted for Model 1 like 
SN magnitude evolution, the Qm parameter is biased when- 
ever |Am evo |i.7| > 0.04. We concluded in the previous section 
that only Am evo |i. 7 < -0.18 and Am evo t i. 7 > 0.14 can be de- 
tected, cf. Table [3] Our ignorance of the evolution thus leaves 
us with a dangerous zone where an SN magnitude evolution 
undetectably biases the results. A full study of all the {/?,} for 
models 1 and 2 yields the results of Table [5] which gives the 
validity zones for Qm, wq and w a , for the four models. 



Table 5. The Am™ 1 1 -7 -intervall in which the fit parameters are 



reconstructed without biases. 









Validity Zone 




Model 








w a 


1 


[ 


-0.04,0.04] 


[-0.08,1.18 ] 


[-0.28,0.15] 


2.1 


[ 


-0.04,0.04] 


[-0.08,0.13 ] 


[ -0.27,0.27 ] 


2.2 


[ 


-0.03,0.05 ] 


[ -0.10,0.37 ] 


[-0.39,0.09 ] 


2.3 


[ 


-0.06,0.09 ] 


[-0.17,0.09] 


[ -0.06, 0.06 ] 



We see that the parameter Qm has the highest chances of being 
biased for all models. Its validity zones are the smallest in com- 
parison to the wq and w a zones. All parameters {Q M , wq, w a }, 
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however, share the general tendency: the higher the power a of 
models 2, the larger the validity zone. A particular feature of 
model 1 is that wq is reconstructed without bias for all positive 



fiducial Am e 



_j < 1.18. 



4.5.3. Danger 

We can merge the obtained detectability and bias-risk limits, 
taken from Tables[3]and[5] respectively, to estimate the danger 
of a neglected magnitude evolution governed by the models of 
Eqs. ( [Tol l and ( fTTT ). We limit our discussion to the parameter 
Qm, because it is the parameter with the highest risk of unde- 
tected biased reconstruction for all models. Figure [5] shows the 
68% confidence level detectability and bias limits for the four 
models obtained from Tabels[3]and[5]in a magnitude-redshift- 
diagram. Dotted lines are the detectability limits obtained by 
the x 2 -test, dashed lines are the detectability limits obtained by 
the most constraining consistency check (underlined values in 
table O, and full lines are the bias limits on Qm taken from 
Table H 

The cross-shaded zones between the detectability and the bias 
limits are the dangerous regions of the diagram, where the 
fiducial evolutionary effect passes the analysis undetected, and 
where biases are introduced at the same time on the parame- 
ter Qm- We may take the area S of the cross-shaded region, in 
the redshift range < z < 1.7, as a measure of the danger of 
the model. The higher the value of S , the greater the risks of 
misinterpretation. We obtain the values of Table[6] 



The analysis of the 95% confidence level naturally yields 
higher absolute values for S, but the same conclusions hold. 
In particular, we obtain the same tendencies within the models, 
and 5 2 > S i for all models. 

The absolute value of S also depends on the specifications of 
the data survey, i.e. in our case on errors on R and A, and on the 
statistics of the SNe in the redshift range. For example, when 
we neglect the assumed systematic error o" sys = 0.02 and de- 
crease the intrinsic magnitude error on each SN to 0. 10, we ob- 
tain bias limits that are tighter than those of Table[5]roughly by 
a factor 2. Also the dectectability limits from Table [3]decrease 
considerably from 30% to 50%, and the^- 2 test becomes more 
and more effective the smaller the intrinsic errors on magni- 
tudes. The absolute values of S also depend on the priors of the 
fit. They may change when the flatness constraint is dropped 
or when w x = const, is imposed on the fit. In the latter case, 
we find a decrease in the detectability ranges and also the bias 
limits. In particular, when im posing w x _ const, o ne obtains re- 
sults consistent with those of ISarkar et alJ (I2008al) . i.e. Qm and 
w x are biased at l<x for evolutionary effects Am evo |i.7 = 0.02 
even in the Model 1 case. 

The measure of danger S we introduce in this section therefore 
should not be taken as an absolute measure, because its actual 
values depend on the specifications of the cosmological probes 
applied in the fit, and on priors. However, we checked that vari- 
ations in the data model or of the priors do not change the con- 
clusion that the late-epoch model Am evo is the most dangerous 
one. 



Table 6. Danger S of the four models as described in the text. 



Model 


sr 


SY 


5 


1 


0.09 


0.14 


0.23 


2.1 


0.17 


0.24 


0.41 


2.2 


0.07 


0.10 


0.17 


2.3 


0.02 


0.03 


0.05 



* Size of the dangerous zone in the upper (lower) part of the Hubble 
diagram only, cf. Fig. [5] 



These values validate the tendency of higher power a ^ 1 evo- 
lutionary effects being less dangerous than a % \ models. The 
late-epoch model 2.1 is more likely to undetectably bias the fit 
results than the early-epoch model 2.3. The danger of the lin- 
ear model 2.2 and the logarithmic model 1 lie in between the 
two extremes, having approximately the same values. Figure|5] 
also shows the values S\ (S 2) of the area of the upper (lower) 
dangerous region seperately. One finds 5 2 > S 1 for all models, 
which again means that a negative magnitude drift Am evo < 
carries the greater risk of undetectably biasing the fit results 
than a positive one. 



We al so e xtracted the danger of th e models used by Linder 
d2006l) and lFerramacho et alJ d2009l) : S = 0.36, which turn out 
to be equivalent, and very close to our Model 2.1. This agrees 
with these parameterizations being well described by a model 
2 with a as 0.6. 



5. Discussion 

We studied four different one-parametric models of SN magni- 
tude evolution on cosmic time scales and obtained constraints 
on its parameters by combined fits on the actual real data 
coming from Supernova surveys, observations of the cosmic 
microwave background, and baryonic acoustic oscillation. 
We found by a minimization of the x 2 that data prefer a 
magnitude evolution of SNe type la such that high-redshift 
supernovae are brighter than would be expected in a standard 
{Q,,w x } cosmos: Am™ < 0. Data are, however, consistent 
with nonevolving magnitudes at the lcr level except for special 
cases. The special cases we found are a fit with the early-epoch 
model 2.3 when dark energy is allowed to be dynamical, 
w x - wo + w- d z/ (1 + z), and a fit with a logarithmic magnitude 
evolution Model 1 when a constant equation of state of dark 
energy w x = const, is assumed in the fit. A comparison of the 
SN magnitude distribution with the CPL best-fit magnitudes 
in a Hubble-diagram also indicate more luminous SN events 
at higher redshift. Our results are, however, of limited strength 
because we do not use the full CMB and BAO data but the 
reduced variables R and A. We also did not include constraints 
coming from other cosmological probes. 
We simulated a future data scenario consisting of 2000 SN 
events out to redshift z = 1.7 and forecasted small errors on 
the CMB shift parameter R and the BAO variable A. In the 
simulation of the SN magnitudes, a redshift dependence of 
the magnitudes according to the four models was allowed for. 
Then in the fit we neglected the possibility of such an evolution 
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and studied the fit results with respect to detectability of this 
wrong model assumption for a wide range of fiducial models. 
We quantified the range of values of the fiducial evolution 
parameters /3 where the wrong assumption is not detectable. 
We found that the linear model 2.2 and the early-epoch model 
2.3 are easier to detect and reject from the fit than is the 
late-epoch model 2.1 and the logarithmic model 1. 

Including the possible biasing of the fitted parameters when 
neglecting magnitude evolution in the discussion, we were able 
to determine the danger of the various fiducial evolutionary 
models; that is to say, we determined the exact parameter 
zones where the wrong model assumption of nonevolving 
intrinsic SN magnitude not only is not detectable, but also 
introduces biases on the fitted cosmological parameters {/?,}. 
The parameter of overall mass density D. M turned out to 
carry the highest risk of biased reconstruction. We found 
that, whereas the dangerous zone is nearly negligible for the 
early-epoch model 2.3 (Am evo ~ z 2 ), it becomes significant in 
the logarithmic model 1 (Am evo ~ log(l + z)) and the linear 
model 2.2 (Am™ ~ z), and dangerous for the late epoch model 
2.1 (Am evo ~ This becomes apparent in Fig. [5] where 
we simultaneously plotted the detectability and bias limits of 
the four fiducial models in a Hubble diagram. A comparison 
of Figs. [5] Q] and |2] shows that the dangerous zone of Model 
2.1 covers nearly the whole area of magnitude dispersion 
of actual SN data in the magnitude-redshift-diagram. Also 
negative magnitude drifts (Am evo < 0) carry the greater risk of 
undetectably biasing the fit results than magnitude evolution 
effects with Am evo > 0, which are already excluded at 1 <x by 
present data for dynamical dark energy. 

In conclusion, special care should be accorded to effects yield- 
ing late-time evolution of SN la magnitudes with a negative 
magnitude drift. Early-time evolution is less severe. To avoid 
any bias in determining the cosmological parameters, it is 
preferable to include a new parameter to describe possible SN 
la magnitude evolution in any combined analysis. From our re- 
sults, it appears that the models chosen to describe such an ef- 
fect should favor late-epoch evolution, as is done by our Model 
2.1, Am evo ~ VE 
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Table 3. The /Hntervall in which an evolutionary effect on SN magnitudes would pass undetected by statistical and physical 
detectability criteria. 



Criterion 


Mc 

68% CL 


del 1 

95% CL 


Mod 

68% CL 


el 2.1 

95% CL 


X 1 < Xcit, 

Consistency 

of: 

< " v 

w a *~* w a r 


[ -0.24 

Zcut = 0.95 

[ -0.20,0.18 ] 
[-0.39,0.33 ] 
[ (...) ,0.41 ] 


,0.26] 

Zc„, = 1.20 

[ -0.18,0.15 ] 
[-0.17,0.13] 
[ -0.22,0.15 ] 


[ -0.28 

Zcut = 0.95 

[ (-0.45), 0.38 ] 
(...) 
(...) 


0.32] 

Zcut = 1-20 

[ -0.43,0.28 ] 
[-0.40,0.23 ] 
[ (...) ,0.26] 


[ -0.2S 

Zcut = 0.95 

[ -0.23,0.20 ] 
[(-0.6), 0.35] 
[ (...) ,0.42] 


,0.28 ] 

Zcut = 1-20 

[ -0.20,0.17 ] 
[-0.19,0.15] 
[-0.25,0.17] 


[ -0.33 

Zcut = 0.95 

[(-0.6), 0.37 ] 
[ (...) ,0.44] 
(...) 


,0.33 ] 

Zcut = 1-20 

[ -0.45,0.32 ] 
[ -0.44,0.26 ] 
[ (...) ,0.30] 



Criterion 


Mo 

68%. CL 


lei 2.2 

95% CL 


Mod 

68% CL 


el 2.3 

95% CL 


X L <xiut 

Consistency 
of: 

Wq *-» nfc r 
w a *~* w a r 


[ -0.12 

Zcut = 0.95 

[ -0.10,0.11 ] 
[-0.20,0.26 ] 
[(-0.6), 0.36] 


,0.14] 

Zcut = 1.20 

[ -0.09,0.08 ] 
[-0.09,0.07 ] 
[-0.12,0.08] 


[ -0.14, 

Zcut = 0.95 

[ -0.21,0.23 ] 
[ -0.38,0.43 ] 
[ ( .) ,(0.7)] 


0.16] 

Zcut = 1.20 

[ -0.21 ,0.14 ] 
[-0.21,0.13 1 
[(-0.3), 0.15] 


[ -0.05 

Zcut = 0.95 

[ -0.04,0.06 ] 
[-0.09,0.20 ] 
[ -0.24,0.28 ] 


,0.05 ] 

Zcut = 120 

[-0.04,0.04] 
[-0.04,0.04 1 
[-0.07,0.07 J 


[ -0.06 

Zcut = 0.95 

[ -0.08,0.12 ] 
[ -0.16,0.28 ] 
[(-0.4), (0.4)] 


,0.06] 

Zcut = 1.20 

[ -0.08,0.07 ] 
[ -0.09,0.07 ] 
[-0.15,0.09] 



Values in round brackets () lie outside the scan range, and were obtained by extrapolation. (...) means that no limits fell within the scan range 
and no meaningful extrapolation was possible. 





Fig. 5. 68% confidence level detectability and bias limits, and the dangerous zones in between (cross-shaded region). 



